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FOREWORD 


This  report  describes  the  research  effort  of  the  Data 
Systems  Division  of  Litton  industries,  Inc.  under  Supple¬ 
mental  Agreement  2  to  Contract  DAAG05-70-C-0328, 
with  the  U  S.  Army,  Frankford  Arsenal.  The  objective 
was  to  provide  additional  analytic  and  simulation  effort  in 
support  of  the  parametric  analysis  of  predicted  fire  air 
defense  weapon  systems 

The  report  is  presented  in  three  volumes.  Volume  I, 
Analysis,  by  Herbert  K.  Weiss,  reports  the  analysis  effort 
and  the  simulation  results.  Volume  II,  Simulation  Model, 
by  Martin  P.  Ginsberg,  describes  the  Litton  Air  Defense 


Simulation,  designed  by  Mr.  Ginsberg.  The  results  of  the 
simulation  are  included  in  Volume  I.  Volume  III,  Effective¬ 
ness,  by  Herbert  K.  Weiss,  reports  on  methods  of  evaluating 
overall  system  effectiveness. 

In  the  present  report,  frequent  reference  is  made  to  the 
Final  Report  on  the  original  contract.  The  previous  report, 
titled  Final  Report,  A  Parametric  Study  of  Advanced 
Forward  Area  Air  Defense  Weapon  System  (AFAADS) 
(two  volumes),  dated  2  October  1 970,  Revised  Edition  1971, 
is  referred  to  throughout  this  report  as  the  "AFAADS-I 
Report." 
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SECTION  1 
INTRODUCTION 


The  AFAADS  gun  model  provides  the  capability  of 
estimating  the  performance  of  proposed  (or  existing) 
predicted  fire  weapon  systems  against  attack  aircraft  in 
one-to-one  combat  situations. 

1.1  CAPABILITY 

The  model  is  not  specific  to  one  gun  or  aircraft 
system  but  has  the  ability  to  model  a  variety  of  sys¬ 
tems.  Each  major  subsystem  of  the  weapon  is  repre¬ 
sented.  The  particular  attributes  or  method  of  opera¬ 
tion  of  each  subsystem  for  a  computer  experiment  is 
determined  by  parameter  values  in  the  model  data  set. 
In  its  present  configuration  the  model  includes  the 
facilities  to: 

a.  Model  realistic  aircraft  paths. 

b.  Model  sensor  behavior  including  the  effects  of : 

( 1 )  Sensor  lag. 

(2)  Regenerative  tracking. 

(3)  Aircraft  image  in  the  tracker  plane. 

(4)  Sensor  noise  due  to  scintillation. 

(5)  Radar  wave  length. 

(6)  Sampling  rate. 

c.  Model  the  data  smoother  for: 

(1)  Position. 

(2)  Velocity. 

(3)  Acceleration. 

d.  Model  the  following  prediction  methods: 

(1)  Linear. 

(2)  Linear  plus  energy  corrections. 

(3)  Linear  plus  acceleration  corrections. 

(4)  Linear  plus  thresholded  acceleration  correc¬ 
tions. 

(5)  Defense  of  a  known  point  combined  with  any 
of  the  above. 

(6)  Simple  polar. 

e.  Model  the  ballistics  of  any  projectile  with  trajec¬ 
tory  characteristics  suitable  to  representation  as  a 
Siacci  tabic  of  no  more  than  300  velocity  points. 

f.  Automatically  use  the  selected  ballistics  and  pre¬ 
diction  algorithm  in  the  time  of  flight,  shell 
position,  and  effectiveness  computations. 

g.  Model  gun  motion  including  the  effects  of: 

( 1 )  Servo  lag. 

(2)  Regenerative  tracking. 


h.  Include  the  effects  of  up  to  10  different  disper¬ 
sion  patterns  in  the  kill  probability. 

i.  Compute  the  single-shot  kill  proba  <ty  for 
equally  spaced  time  points  over  the  aircrmi  path. 

j.  Compute  the  burst-kill  probability  for  each  burst 
over  the  aircraft  path  for  up  to  10  fire  doctrines 
on  each  execution. 

k.  For  each  fire  doctrine  and  dispersion  combina¬ 
tion  compute  the  probability  of  a  kill. 

l.  Produce  detail  reports  of  the  performance  of  any 
subsystem. 

m.  Produce  summary  reports  of  each  execution. 

n.  Automatically  replicate  and  produce  summary 
statistics  across  replicates  of: 

(1)  The  distribution  function  of  the  miss-dis¬ 
tances  observed  over  the  course. 

(2)  The  distribution  function  of  the  single-shot 
kill  probabilities  over  the  course  for  each 
dispersion  pattern  input. 

(3)  The  kill-sec  (integral  or  single-shot  kill  proba¬ 
bility)  over  the  course  for  each  dispersion 
pattern  input. 

(4)  The  distribution  function  of  burst  kill  proba¬ 
bilities  for  each  dispersion  pattern,  fire  doc¬ 
trine  combination  input. 

(5)  The  kill  probability  over  the  course  for  each 
fire  doctrine  and  dispersion  set  averaged  over 
the  set  of  replicates. 

1.2  EXPERIENCE 

Figure  I- 1  is  a  top  level  block  diagram  of  the  model 
organization.  The  model  was  intentionally  structured 
for  use  as  a  study  tool.  The  area  enclosed  in  dashed 
lines  is  where  replication  and  statistics  gathering  over 
samples  using  the  same  data  values  and  different  ran¬ 
dom  number  streams  are  accomplished.  The  outer  loop 
which  encloses  this  area  is  used  for  the  selective  re¬ 
placement  of  parameter  values.  This  structure  enables 
the  completion  of  an  experiment  in  one  computer  run. 
For  example,  if  it  is  desired  to  construct  a  footprint, 
the  inner  loop  would  compute  statistics  over  replicas 
for  one  path  and  the  outer  loop  would  be  used  to 
change  the  minimum  crossing  distance  of  the  path. 

The  AFAADS  model  has  been  used  to  conduct  com¬ 
puter  experiments  to  measure  the  variability  of  system 
performance  due  to  changes  in  aircraft  performance, 
system  implementation,  and  system  operating  policy. 
Specific  experiments  which  have  been  conducted  and 
reported  include  effects  of : 
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•  Varying  aircraft  path  including: 

Varying  velocity  on  linear  paths. 

Varying  crossing  range  on  linear  paths. 

Introducing  jinking  on  linear  paths. 

Introducing  'break-away*  maneuvers  on  linear 
and  jinking  paths. 

•  Varying  the  sensor  subsystem  with  respect  to: 

The  wave  length  of  radar. 

Using  manual  tracking. 

Using  rate-aided  manual  tracking. 

Introducing  regeneration 
Introducing  angular  bias. 

Interruption  of  sensor  data. 

•  Varying  the  smoothing  algorithms  with  respect 
to: 

Smoothing  period. 

Number  of  points. 

•  Varying  the  prediction  algorithm  including: 
Linear  prediction 

Quadratic  prediction. 

Energy  correction. 

Defense  of  a  known  point. 

•  Varying  the  projrctile-to-projectile  position  devi¬ 
ation  due  to: 

Angular  dispersion. 

Velocity  dispersion. 

•  Varying  the  ballistics  including: 

Muzzle  velocity. 

Caliber. 

•  Varying  the  fire  doctrine  including: 

Continuous  vs.  interrupted  fire. 

Rate  of  fire. 

Open  fire  range 
Burst  duration. 

In  addition  to  the  studies  listed  above  the  model  has 
also  been  used  to  generate  •footprints'  for  20,  25,  37 
and  50  mm  gun  systems  against  aircraft  executing  both 
unaccelerated  and  jinking  fly-by  courses. 

The  simulation  model  is  programmed  in  FORTRAN 
IV  and  is  currently  executing  on  IBM  370/165  and 
CDC  6400  computer  systems.  Execution  times  are 
highly  dependent  on  the  particular  experiment  being 
conducted,  however,  the  model  is  fairly  fast.  For  exam¬ 


ple,  to  execute  a  fly-by  path  from  acquisition  at  3600 
meters  to  crossover  for  10  fire  doctrines  and  S  disper¬ 
sion  patterns  requires  less  than  I  CPU  second  per 
replicate  on  the  IBM  370/I6S. 

1.3  MODEL  ELEMENTS 

The  exercise  portion  of  the  model  is  structured  to 
maintain  the  'real  world’  information  pattern  of  the 
gun  system  being  simulated.  This  is  accomplished  by 
restricting  the  access  of  each  simulated  gun  component 
to  only  that  information  which  would  be  available  in 
the  'real  world.’ 

The  manner  in  which  information  is  acquired,  gen¬ 
erated.  and  transmitted  through  the  system  is  pre¬ 
sented  in  Figure  1-2.  The  primary  information  flow  is 
as  follows:  the  aircraft  position  as  a  function  of  time 
enters  the  system  via  the  Data  Acquisition  Event, 
which  adds  the  effects  of  the  sensor  system  and  fills  the 
Sensed  Position  File;  the  Prediction  Event  uses  the 
sensed  data,  the  time  of  flight  module,  and  the  predic¬ 
tion  algorithm  module  to  fill  the  Predicted  Position 
File  with  aim  points;  the  Gun  Motion  and  Fire  Event 
accesses  the  Predicted  Position  File  and  fills  the  Shell 
Position  File  as  if  the  gun  fired  continuously;  the 
Report  Generator  Event  accesses  the  Shell  Position 
File  and  the  Aircraft  Position  Module,  applies  the  fire 
doctrines  called  for,  and  produces  reports.  The  preced¬ 
ing  information  flow  pattern  demonstrates  an  impor¬ 
tant  feature  of  the  model.  Note  that,  as  in  the  ‘real 
world’,  any  action  which  ihe  gun  system  lakes  (or  can 
take)  is  ultimately  based  on  information  acquired 
through  its  sensing  system 

Table  I- 1  lists  the  major  elements  of  the  computer 
model.  Each  is  described  in  greater  detail  in  the  fol¬ 
lowing  paragraphs  and  specific  formulation  and  equa¬ 
tions  for  the  elements  arc  presented  in  the  main  body 
of  the  text. 

1.3.1  Aircraft 

The  aircraft  is  described  by  its  flighi  path,  its  size, 
and  its  shape.  The  flight  path  is  specified  as  an  initial 
position,  velocity,  and  heading,  and  a  series  of  path 
segments  or  legs.  The  legs  are  described  by  specifying 
their  horizontal  and  vertical  projections.  Horizontal 
projections  may  be  segments  of  lines  (accelerated  or 
unaccelerated),  segments  of  circles,  or  segments  of 
spirals.  Vertical  projections  may  be  segments  of  lines 
or  segments  of  parabolas.  The  bank  angle  of  the 
aircraft  is  determined  by  the  radial  acceleration  of  the 
segment  in  question. 

Aircraft  size  and  shape  are  considered  in  the  data 
acquisition  phase  and  in  the  vulnerability  phase.  In  the 
data  acquisition  phase  the  aircraft  is  represented  as 
three  perpendicular  vectors:  The  fuselage  length  vec¬ 
tor,  along  the  line  of  flight;  the  wing  length  vector, 
perpendicular  to  the  fuselage  length  vector  and  rotated 
from  the  horizontal  plane  by  an  angle  equal  to  the 


Table  1-1.  Major  Model  Elements 
i  AIRCRAFT 

Might  Pith 
Sire  and  Shape 

II  WFAPON  SYSTEM 

Sensor 

Smoother 

Predictor 

Controller 

Gun 

Projectiles 

III  RhPORTCI  NKRATOR 

Aircraft  Vulnerability 
I  ire  Doctrine 

Measures  of  Performance  for: 

Sensor  Smoother 
Predictor 
Controller 
Weapon  System 
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bank  angle;  and  the  aircraft  height  vector,  in  a  direc¬ 
tion  perpendicular  to  the  first  two  vectors  The  lengths 
of  these  vectors  are  combined  with  the  aircraft-sensor 
geometry  in  determining  the  standard  deviations  of 
the  sensor  error  in  azimuth  and  elevation. 

Aircraft  size  and  shape  for  the  vulnerability  phase 
are  entered  separately.  The  aircraft  is  represented  as 
two  perpendicular  ellipsoids:  one  for  the  wing  and  one 
for  the  fuselage.  Each  ellipsoid  is  specified  by  its  three 
semi-major  axes:  length,  width,  and  height.  Vulnerabil¬ 
ity  of  the  aircraft  then  is  determined  by  combining  the 
vulnerability  of  the  fuselage,  the  wing,  and  the  fuse¬ 
lage-wing  intersection 

1.3.2  Weapon  Syatam 

The  weapon  system  sensor  model  uses  the  actual 
aircraft  position  as  a  departure  point  and  adjusts  the 
values  to  reflect  aircraft-sensor  geometry,  tracker  drive 
elements,  and  stochastic  error  due  to  ‘glint’. 

The  actual  position  is  obtained  in  cartesian  coordi¬ 
nates  and  transformed  to  polar  coordinates.  The  posi¬ 
tion  and  velocity  of  the  aircraft  art  used  to  determine 
the  size  of  the  aircraft  in  the  sight  plane.  These  sizes 
are  combined  with  azimuth  and  elevation  rates  to 
determine  the  size  of  the  error  distributions  in  azimuth 
and  elevation,  and  the  current  amount  of  error  correla¬ 
tion  across  time.  The  history  of  tracker  motion  is  used 
lo  determine  the  amount  of  tracker  lag  and  past  sensed 


values  are  used  to  determine  how  much  regeneration  is 
present  (if  specified).  Random  numbers  are  used  to 
determine  current  values  for  azimuth,  elevation,  and 
range  errors  due  to  glint.  The  adjusted  values  are 
transformed  to  cartesian  coordinates  and  stored  as 
sensed  position. 

The  weapon  system  data  smoother  operates  on 
sensed  data.  The  particular  filters  to  use  are  determined 
by  input  data.  Three  filters  arc  required:  Position, 
velocity  and  acceleration.  The  sensed  position  file 
maintains  the  ten  most  recent  samples  at  all  times. 
Thus,  any  filter  which  can  be  represented  by  ten  or 
fewer  coefficients  to  operate  on  equally  spaced  data 
values  can  be  used. 

The  weapon  system  predictor  operates  on  smoothed 
data.  There  are  several  algorithms  available  with  the 
specific  one  used  on  z  given  experiment  determined  by 
the  input  data.  Algorithms  currently  available  are: 

a.  LINEAR.  Straight  linear  extrapolation  of  current 
smoothed  position  based  on  current  smoothed 
velocity  and  time  of  flight. 

b.  QUADRATIC.  Quadratic  extrapolation  of  cur¬ 
rent  smoothed  position  based  on  current 
smoothed  position,  velocity  and  radial  accelera¬ 
tion. 

c  QUADRATIC  THRESHOLD.  Uses  linear  ex¬ 
trapolation  unless  the  radial  acceleration  exceeds 
a  specified  threshold.  Automatically  adds  an  ac¬ 
celeration  correction  to  the  linear  prediction 
whenever  the  threshold  is  exceeded. 

d.  SIMPLE  POLAR.  Linear  extrapolation  of  air¬ 
craft  position  in  polar  coordinates  based  on  mea¬ 
sured  position  and  rates. 

e  DEFENSE  OF  A  KNOWN  POINT.  Assumes  air¬ 
craft  will  turn  and  attack  a  known  point  if 
current  heading  and  velocity  allow  such  a  ma¬ 
neuver.  Used  in  conjunction  with  one  of  the 
above  algorithms. 

f.  CONSTANT  ENERGY.  In  conjunction  with  the 
above  algorithm  a  velocity  correction  based  on 
velocity  in  the  vertical  plane  can  be  used.  It  is 
based  on  a  thresholding  mechanism  similar  to 
the  automatic  linear-quadratic  switch. 

Aim  errors  due  to  the  gun  subsystem  are  of  two 
types:  lag  between  aim  order  and  actual  arrival  at  the 
aim  point,  and  bias  due  to  imperfections  in  the  system 
alignment  and/or  adjustment.  Lag  is  computed  using  a 
history  of  the  commands  given  to  the  gun  system 
Regeneration  may  be  specified  and  when  present  is 
also  based  on  the  command  history.  Provision  has  been 
made  for  the  specification  of  systematic  errors  in 
azimuth  and/or  elevation. 
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Projectile  trajectories  are  computed  using  a  set  of 
input  tables.  Any  family  of  trajectories  which  can  be 
specified  using  tables  (of  less  than  500  points  each) 
expressing  time  of  flight  and  range  as  a  function  of 
remaining  velocity  may  be  used.  Provision  has  been 
made  for  the  specifics  aon  of  a  bias  of  the  muzzle 
velocity. 

1.3.3  Report  Generator 

The  report  generator  portion  of  the  model  is  much 
more  than  an  output  routine  For  economy  of  com¬ 
puter  tinr  certain  computations  are  made  in  the 
report  ge  rator,  specifically,  the  effects  of  both  round- 
to-roum  ispersion  and  fire  policy,  as  well  as  aircraft 
vulnerai.iity.  are  computed  in  the  report  generator. 
This  arrangement  allows  the  generation  of  reports  of 
up  to  10  fire  policies  and  10  sets  of  dispersion  parame¬ 
ters  using  the  results  of  one  mode!  execution. 

During  the  exercise  portion  of  the  model,  the  gun  is 
‘fired’  continuously,  in  the  Report  Generator,  fire  doc¬ 
trines  are  used  to  determine  when  the  gun  system 
would  fire,  how  many  rounds  and  the  rate  of  fire 
Current  implementation  requires  that  the  following 
data  he  specified: 

a.  Amount  of  ammunition  on  the  mount 

b.  Number  of  rounds  per  burst 

c.  Open  fire  time 

d.  Duration  of  a  burst 

e  Dead  time  (time  between  bursts)  and 

f.  Reload  time. 

The  appropriate  miss  di;  lances  and  the  round-to- 
round  dtsoerston  parameters  are  combined  with  the 
aircraft  size  in  the  VULNERABILITY  module  tc  pro¬ 
duce  kill  probabilities. 


Measures  of  performance  may  be  generated  for 
subsystem  or  the  gun  system.  They  may  be  at  the  detail 
level  for  each  replicate  or  at  the  summary  level  over  a 
set  of  replicates. 

1.4  STRUCTURE  OF  THE  REPORT  BODY 

In  the  following  sections,  the  model  capability  is 
described  in  greater  detail.  To  facilitate  the  discussion, 
ihe  model  is  viewed  in  three  sets  of  functionally  re¬ 
lated  program  modules 

The  first  set  of  modules  is  used  to  control  the  simula¬ 
tion  process.  This  set  is  composed  of  the  main  program 
and  one  subroutine  which  maintains  a  schedule  of 
events  and  provides  for  central  time  accounting. 

The  second  set  of  modules  is  composed  of  the  Gun 
System  Functional  Event  Routines.  This  set  of  routines 
models  the  various  components  of  the  gun  system 
under  study.  There  are  five  event  routines: 

a  Data  Acquisition  Event. 

b.  Sensor  Failure  Event. 

c  Prediction  Event 

d  Gun  Motion  and  Firing  Event. 

e  Report  Generation  Event. 

The  third  set  of  modules  is  the  utility  package.  These 
modules  are  used  by  the  event  routines  or  each  other 
and  may  not  be  invoked  by  the  system  scheduler. 
Included  in  this  set  are: 

a.  Aircraft  positioner. 

b.  Time  of  flight. 

c.  Prediction  algorithm. 

d  Random  number  generator 

c.  Coordinate  transformation  package 

f.  Fire  Doctrine  Module. 

g.  Vulnerability  Computation  Module. 
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SECTION  2 

SIMULATION  CONTROL 


The  simulation  model  represents  the  combat  situa¬ 
tion  under  study  as  a  series  of  events  centered  around 
the  antiaircraft  gun  system.  Event  sequencing  and 
initialization  are  performed  by  the  master  control 
program  and  the  clock  module.  The  initialization 
phase  consists  of  reading  input  parameters,  perform¬ 
ing  and  storing  the  results  of  once-only  computations, 
and  scheduling  the  first  (a  sensor  data  acquisition 
event)  and  last  (a  report  event)  evenu  to  be  processed. 
Initialization  tasks  are  accomplished  by  special  mod¬ 
ules  which  are  called  by  the  master  control  program 
prior  to  entering  the  event  sequencing  mode. 

The  logic  of  the  control  process  is  presented  in 
Figure  2-1.  The  particular  order  of  evenu.  which  will 
occur,  is  determined  by  the  gun  system  and  the  aircraft 
path  under  study.  There  are.  however,  some  principles 
worth  noting:  for  example,  the  first  event  to  occur  will 
always  be  Event  Type  I.  and  it  will  continue  to  occur 
at  its  specified  period  (input  during  initialization 
phase)  until  the  end  simulation  event.  With  the  excep¬ 
tions  of  the  end  simulation  event  and  the  sensor  inter¬ 


ruption  evenu,  no  other  evenu  will  occur  unless  they 
are  explicitly  scheduled  by  Event  Type  I  or  some  other 
module,  either  directly  or  indirectly  scheduled  by  an 
occurrence  of  Event  Type  1 .  The  occurrence  of  Event 
Type  2  indicates  the  end  of  the  current  sequence  of 
simulation  evenu. 

One  of  two  actions  are  initiated  after  the  occurrence 
of  Event  Type  2.  If  there  are  more  replications  to  be 
performed;  the  event  file  is  cleared,  all  Dags  are  turned 
off,  the  initial  Event  Type  I  and  end  simulation  evenu 
are  scheduled,  and  the  simulation  process  continues.  If 
no  more  replications  are  required,  the  summary  statis¬ 
tics  over  replications  are  printed,  and  the  model  enters 
the  selective  input  phase.  During  the  selective  input 
phase  any  data  set  may  be  replaced  and  the  model 
restarted.  The  selective  input  phase  is  very  useful  for 
making  sensitivity  studies  since  a  series  of  exercises 
can  be  made  specifying  the  full  data  set  once  and  only 
the  data  which  changes  from  case  to  case  thereafter. 
This  enables  the  completion  of  an  extended  computer 
experiment  in  one  computer  run. 


2-1 


MORE 

REPLICATES 


REPORT 

GENERATION 


MORE 

EXERCISES 


SELECTIVE 

INPUT 

INITIALIZATION 


F'flure  2-1.  Logic  of  the  Control  Process 
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SECTION  3 
EVENTS 


3.1  SENSOR  DATA  ACQUISITION  EVENT 

In  the  real  world  all  gun  system  functions,  which 
required  a  knowledge  of  aircraft  position  or  rate  pa¬ 
rameters,  acquire  the  information  either  directly  or 
indirectly  from  the  gun  sensor.  In  the  computer  model 
the  sensor  function  is  modeled  by  the  Sensor  Data 
Acquisition  Event.  On  each  entry,  this  event  routine 
stores  current  values  of  sensed  aircraft  position  in  the 
SENSED  DATA  File.  The  gun-system’s  data  smoother 
is  also  modeled  in  this  routine  and  its  outputs, 
smoothed  position  and  rate  values,  are  also  placed  in 
the  data  file. 

In  the  following  discussion,  upper  case  letters  will  be 
used  for  true  values,  primed  upper  case  letters  for 
sensed  values,  underscored  upper  case  letters  for 
smoothed  values  of  position,  and  dot  notation  to  indi¬ 
cate  rates  and  accelerations.  The  time  to  which  values 
pertain  will  be  indicated  by  subscript  expressed  relative 
to  current  time.  Thus,  if  A  is  a  position  coordinate: 


A, 

m 

Current  true  value. 

A'. 

m 

Current  sensed  value 

A. 

m 

Current  smoothed  value. 

A. 

m 

Current  smoothed  rate  value. 

A' t 

m 

Sensed  value  at  current  time  minus  At. 

The  six  dimensional  transformations  from  cartesian 
to  polar,  and  from  polar  to  cartesian  coordinates,  will 
be  represented  by  the  operators;  ‘CTOP’  and  ‘PTOC.' 
Thus: 


(D,  6.0,  b,0,<p)  *  CTOP  (X.  Y.  Z.  X.  Y.  Z)  (3. 1 ) 
and 

(X.  Y,  Z.  X,  Y,  Z)  -  PTOC  ID.  6,  0,  6. 6 , 0\  (3.2) 


On  each  entry  to  the  event  routine,  the  sequence  of 
computations  is  as  follows.  The  true  aircraft  position 
and  rates  are  obtained  by  invoking  the  aircraft  posi¬ 
tion  module.  These  values  which  are  returned  in  carte¬ 
sian  coordinates  are  transformed  to  polar  coordinates, 
which  is  the  natural  coordinate  system  of  the  sensor 
systems.  From  the  sensor-aircraft  geometry,  the  size  of 
the  aircraft  as  ‘seen’  in  azimuth,  elevation,  and  range 
are  computed  and  saved.  Based  on  the  aircraft  size  and 
rates  in  the  sight  plane,  values  are  computed  for  the 
standard  deviation  of  glint  noise  in  azimuth  and  eleva¬ 
tion  and  for  the  sample-to-sample  correlation.  Current 
aircraft  position  is  estimated  based  on  the  last  values 
of  smoothed  position  and  velocity.  A  check  is  made  to 
determine  the  state  of  the  sensor.  Any  values  which  are 
in  an  interrupted  state  are  replaced  by  the  values 
computed  from  the  smoothed  rates  and  further  adjust¬ 
ments  to  these  values  are  not  made.  Values  in  ar.  up 
state  are  adjusted  for  lag,  regeneration,  and  the  pres¬ 
ence  of  sensor  noise. 

The  resulting  sensed  position  in  polar  coordinates  is 
transformed  to  cartesian  coordinates  and  stored  in  the 
SENSED  DATA  File.  The  most  recent  ’n*  values  stored 
in  the  SENSED  DATA  File  are  used  to  compute 
current  values  for  smoothed  position,  velocity,  and 
acceleration.  These  values  are  also  stored  in  the 
SENSED  DATA  File. 

The  dimensions  of  the  aircraft  in  the  sight  plane  are 
obtained  by  a  rotation  and  projection  onto  the  sight 
plane  of  a  set  of  three  orthogonal  vectors  representing 
the  aircraft  length,  height,  and  width.  First,  a  rotation 
to  orient  the  set  in  space  such  that  the  length  is  parallel 
to  the  aircraft  velocity  vector  and  the  width  is  rotated 
to  an  angle  equal  to  the  bank  angle  is  performed.  Then 
these  rotated  vectors  are  projected  onto  the  sight  plane 
which  is  perpendicular  to  the  range  vector  and  the 
maximum  value  of  azimuth,  and  elevation  projections 
are  determined  and  saved. 


3.1.1  Senior  loo 


Standard  superscript  arrow  notation  (A)  will  be  used 
to  designate  vectors.  Thus,  the  coordinate  transforma¬ 
tions  become: 


The  computation  for  sensor  lag  is  based  on  actual 
aircraft  rates  in  polar  coordinates,  since  these  are  the 
cause  of  the  phenomenon  in  the  Teal  world.’  Lag  is 
computed  in  the  azimuth  and  elevation  coordinates 
using  the  following  operation: 


and 


T-  CTOP  (cj 

(3.3) 

«  • 

L  -  1  p  +-L  P°'£a 

0  K/o  Ka  A 

PTOC  fP) . 

(3.4) 

This  is  the  usual  two  term  servo  lag  approximation 
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with  acceleration  being  estimated  using  two  rate  mea¬ 
surements 

3.1.2  Regenerated  Rates 

The  computation  of  rates  for  regenerative  sensor 
tracking  is  based  on  data  acquired  by  the  sensor.  This 
is  again  in  agreement  with  the  information  flow  in  the 
‘real  world*  situation.  Estimates  of  current  position  are 
computed  from  current  values  of  the  data  smoother 
using  linear  extrapolation: 

xS-*.a+2.aa 

ZJ-Z^+Z.^  A  (3.6) 


where: 

X,  -  is  the  current  value  of  sensor  noise  in  the 
coordinate 

X,t  -  is  the  last  value  of  sensor  noise  in  the 
coordinate 

N  •  is  a  sample  from  a  population  which  is 
normally  distributed  with  mean  zero  and 
variance  one 

0  -  are  coefficients  to  obtain  the  correct  serial 
correlation 

o  •  is  the  standard  deviation  of  the  noise  dis¬ 
tribution  in  the  coordinate. 

For  the  range  coordinate,  a,  fi,  and  o  are  static  over 
a  computer  run.  The  value  of  o  is  an  input  parameter 
and  a  and  p  are  computed  as  follows: 


Ground  range.  S.  and  ground  range  rate.  S,  are 
estimated  as  follows: 

S-(X*2  +  Y*2)1/2  (3.7) 

S»(X*£.A*Y  (3.8) 

Polar  position  and  rate  are  obtained  as  follows: 

P*  -  CTOP  (X*.  V0,  Z*.  X  .A,  Y  i  .4)  (3.9) 

Angular  accelerations  are  computed  from  the  above 
values  by: 


0*»-2.O0*S/S 


(3.10) 


6*  »  -2.0  (D*D*0*  +  S0*2Z*)/D*2  (3.1 1 ) 


0-(i  -  a2) 1/2 


(3.15) 

(3.16) 


where: 


A  is  the  sample  interval 
r  is  the  autocorrelation  of  the  sample. 

For  the  angular  coordinates,  a.  p.  and  o  are  dynami¬ 
cally  computed  using  parameters  which  describe  the 
sensor  and  the  sensor-target  geometry.  They  are  de¬ 
fined  by  the  following  equations: 

|  (i*5oo^)  (|,5rc)  i'/! 

<’"'2  07  )— 7%'wH -bk 


tag) 

(i*,000k)! 


l1* 'Ob?;  2 


Finally,  the  amount  of  regeneration  provided  is 
given  by: 


(3.12) 


Cl 

•B  ,0  t 


.eBA 


.f.0§t>,  (3.18) 

BX 


ROb  +  — 

Kv  Ka 


(3.13) 


3.1.3  8tnsor  Noise 


The  value  for  sensor  noise  in  each  of  the  coordinates 
is  computed  using  the  formula: 


Xj  ■aXj.|  +0<tn 


(3.14) 


where: 


f  -  [I  a2|  1/2 


(3.19) 


f)  is  the  angular  rate  in  the  coordinate  and  is 
composed  of  the  rates  due  to  aircraft  path 
and  ind 

L  is  the  length  of  the  aircraft  in  the 
coordinate 


B  is  the  bandwidth  of  the  servo 


X  is  the  wave  length  of  the  radar 
A  is  the  sample  interval. 

There  are  provisions  in  the  model  for  the  case  when 
there  is  an  interaction  between  lag  and  noise.  When 
this  condition  exists,  the  equation  becomes: 


Xj  *  oXj.,  +  (/to  ♦  *i(L+R)J  N  (3.20) 


where  the  new  term  p(L+R)N  is  used  to  increase  the 
standard  deviation  of  the  noise  sample  by  a  fraction, 
ft,  of  the  existing  lag  (the  sum  of  lag  and  regeneration) 
at  sensing  time. 

3.1.4  Systematic  Errors  (Biss) 

Systematic  errors  in  any  combination  of  range,  azi¬ 
muth,  and/or  elevation  may  be  included.  The  value  of 
b'at  in  each  of  these  coordinates  is  part  of  the  input 
'<4ta  set  and  remain  constant  over  a  computer  run. 

3. 1.5  Sensed  Position 

The  sensed  position  stored  for  a  given  time  is  given 
by  the  equation: 

VPo  +  Lo  +  Ro*No+B  (321) 

where: 

P  indicates  polar  coordinates 
L  is  lag 

R  is  regeneration 
N  is  noise 
B  is  bias 

In  addition,  the  valuet  are  stored  in  the  position  file 
in  cartesian  form.  That  it: 

Cg  ■  PTOC  (Pq)  it  saved  for  use  in  computing 
smoothed  values. 

(3.22) 

3.1.6  Smoothed  Vslues 

The  tracker  computet  and  stores  valuet  for  current 
smoothed  position,  velocity,  and  acceleration.  These 
values  are  computed  using  the  ‘n’  most  recent  values  of 
sensed  position.  In  equation  form: 


c°'h  ''■» 

(3.23) 

n-1 

(3.24) 

-»  "•*  -* 

So-L  'jSa 

J-0 

(3.25) 

3.1.7  Regenerative  Tracking 

The  Sensor  Data  Acquisition  event  hat  the  capability 
of  generating  'tensed’  positions  which  arc  bas.-d  on 
extrapolations  of  smoothed  data  values.  The  event 
automatically  uses  this  method  of  providing  sensed 
data  whenever  data  flow  from  the  sensor  is  inter¬ 
rupted.  The  interruption  can  be  due  to  a  tensor  failure 
in  range,  in  angles,  or  in  both  range  and  angle;  or  ue 
to  terrain  masking.  Sensor  failures  are  scheduled  in  the 
input  phase  and  when  they  occur,  the  appropriate 
parameters  to  in'errupt  the  data  flow  from  the  sensor 
are  set.  Terrain  Masking  is  specified  as  part  of  the 
aircraft  path.  In  the  cate  of  terrain  masking,  ihc 
aircraft  position  routine  returns  an  indicator  to  the 
Data  Acquisition  Event  whenever  the  aircraft  is 
masked,  when  this  indicator  it  on.  the  regenerated 
coordinates  are  used  instead  of  sensed  values. 

Table  III- 1  presents  the  method  used  to  obtain 
’sensed’  values  in  the  presence  of  other  than  normal 
sensor  operation.  The  equations  show  how  the  ’sensed' 
position  is  computed.  Whenever  one  or  more  coordi¬ 
nates  are  still  operable,  lag  and  noise  are  added  to 
them  as  in  normal  operations;  however,  no  adjust¬ 
ments  are  made  to  regenerated  coordinates. 

3.1.8  Event  Scheduling 

The  DATA  ACQUISITION  EVENT  also  performs 
the  initial  scheduling  for  the  PREDICTION  EVENT 
and  the  OUN  FIRING  EVENT.  On  the  first  occur¬ 
rence  of  the  DATA  ACQUISITION  EVENT  with  the 
aircraft  in  a  visible  state,  both  the  PREDICTION  and 
GUN  EVENTS  will  be  scheduled.  The  initial  time  for 
each  of  these  events  will  be  the  time  of  th*  DATA 
ACQUISITION  EVENT  plus  the  system  acquisition 
time,  TACQ,  which  is  an  input  variable. 

3.2  SENSOR  FAILURE  EVENTS 

SENSOR  FAILURE  is  perhaps  an  error  in  naming. 
There  are  four  events  which  are  processed  by  this 
event  routine.  They  are: 

a.  Range  Interruption. 

b.  Range  Resumption. 
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Table  III- 1 .  Regenerated  Coordinate  Computation 


Sute 

Computation 

Range  Failure 

c  ■  PT0C|D*.eo^o,D*.*o,;o| 

Angle  Failure 

Z  ■  PTOC  |D0,fl*,**,D0.$*.  J*| 

Range  and  Angle  Failure 

c  -  5* 

Terrain  Masking 

z-  -e* 

C*  -  |X*.  Y*.  Z*.  X*.  &•,  t*\  • 

12k  ♦  1 *  A*i  x.  5,  ij 

p»  -  |D».flv *•. -  crop t£*i 

2087I-SJ3A 


c.  Angle  Interruption. 

d.  Angle  Resumption. 

Each  it  tchcduled  during  the  initialization  phase  and 
may  occur  periodically  if  the  period  it  lest  than  the 
time  duration  of  the  exercise.  The  total  extent  of  the 
action  taken  by  the  event  when  it  occurs  it  to  set  or 
reset  an  appropriate  variable  to  represent  the  current 
state  of  the  tensor.  The  variable  it  used  in  the  DATA 
ACQUISITION  EVENT  to  determine  the  state  of  the 
sensor  and  perform  appropriate  computations. 

3  3  PREDICTION  EVENT 

The  Prediction  Event  models  that  portion  of  the  real 

!;un  system  in  which  tensor  information  is  extrapo- 
ated  into  the  future  in  order  to  compute  lead  angles 
for  gun  pointing.  T  his  event  routine  can  operate  in  two 
modes,  the  first,  standard  mode,  makes  no  restraining 
assumptions  on  the  aircraft  trajectory.  The  second, 
defense  of  a  known  point,  assumes  that  a  known  point 
is  the  object  of  attack  and  biases  predictions  according 
to  that  assumption. 


The  sequence  of  compulations  which  occur  when  an 
entry  is  made  to  this  event  routine  are  at  follows:  First, 
the  Time  of  Flight  Module  it  invoked  which,  with  the 
aid  of  the  Prediction  Algorithm  Module,  computet  a 
time  of  flight.  The  returned  value  for  time  of  flight  it 
combined  with  current  time  and  the  Prediction  Algo* 
rithm  Module  is  invoked  directly  to  obtain  the  re¬ 


quired  predicted  position.  At  this  point,  if  standard 
prediction  hat  been  specified,  a  branch  is  executed  to 
the  area  which  stores  data  in  the  PREDICTED  POSI¬ 
TION  File.  If  defense  of  a  known  point  has  been 
specified,  the  process  continues. 


As  mentioned  above,  when  defense  of  a  known 
point  is  specified,  predictions  arc  biased.  The  bias  may 
take  two  forms  based  on  either  the  assumption  that  :ne 
aircraft  will  turn  toward  the  defended  point,  or  the 
assumption  that  the  aircraft  will  follow  a  straight  line 
from  its  present  position  to  the  defended  point.  Con¬ 
sideration  of  the  geometry  of  the  situation  reveals  that 
neither  of  these  biases  are  appropriate  under  certain 
conditions.  For  this  reason,  even  when  bias  is  specified, 
a  particular  prediction  may  or  may  not  be  biased.  The 
determining  factor  is  the  position  and  velocity  of  the 
aircraft  relative  to  the  defended  point.  If  the  defended 
point  is  in  front  of  the  aircraft,  a  bias  is  considered.  ‘In 
f ront  of  should  be  understood  as  an  angle  of  less  than 
90  degrees  between  the  aircraft  horizontal  velocity 
vector  and  the  horizontal  projection  of  the  vector  from 
the  aircraft  to  the  defended  point.  Since  the  dot  pro¬ 
duct  of  two  vectors  is  positive,  if  and  only  if  the 
absolute  value  of  the  angle  between  the  vector  involved 
is  less  than  90  degrees,  the  decision  rule  becomes: 

a,  Consider  bias  if:  V  •  R  >  0. 

b.  Do  not  consider  bias  if:  V  •  R  <  0. 


6  ■  V^/ng 


(3.28) 
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Figure  3-1.  Bias  Consideration  Region 

Figure  3-1  displays  the  geometry  of  tue  ueciMon 
rule.  If  the  defended  point  lies  within  the  shaded  area 
bias  is  considered. 

Whenever  the  defended  point  lies  within  the  bias 
consideration  area  of  the  aircraft,  the  geometry  is 
further  examined  to  determine  if  a  bias  can  be  used, 
and  if  so.  which  of  the  two  biases  is  appropriate.  Three 
situations  are  possible;  I )  the  aircraft  heading  is  such 
that  a  projection  of  its  current  velocity  vector  would 
pass  within  a  specified  (during  initialization)  distance 
(R«)  of  the  defended  point.  2)  a  standard  turn  of 
specified  (during  initialization)  acceleration  followed 
by  a  straight  line  flight  path  would  cause  the  aircraft  to 
pass  over  the  target,  of  3)  neither  of  the  above. 

The  first  situation,  namely  the  straight  line  bias,  is 
the  easiest  to  explain  (see  Figure  3-2).  The  heading 
angle  change  0  is  computed  as  follows; 


- 1/2  £t )  (3.26) 


The  closest  approach  to  the  defended  point  is  ap¬ 
proximated  by  the  arc; 


...  _  .|  /RxV 

*“Tan  (iTv 


S  •  Rd  (3.27) 

Linear  bias  is  used  if  S  «  R,. 

The  linear-bias  computations  are  accomplished  in 
the  Time  of  Flight  Module  and  the  Prediction  Algo¬ 
rithm  Module.  If  a  situation  is  such  that  linear  bias  is 
appropriate,  an  indicator  switch  is  set  and  point  to 
point  interpolation  will  be  used  in  the  Prediction 
Algorithm  Module. 

If  the  situation  it  such  that  linear  bias  is  not  appro¬ 
priate.  a  check  is  made  to  see  if  a  turn  followed  oy  a 
straight-line  path  can  bring  the  aircraft  over  the  target. 
This  check  is  made  as  follows:  The  radius  of  turn  is 
computed  using  the  following  equation: 


where: 


V  »  the  aircraft  velocity 
ng  -  ihe  standard  turn  acceleration  (n  is  input) 

The  heading  and  turn  radius  are  used  to  And  the 
center  of  curvature  of  the  tum  and  the  distance  (L) 
from  this  point  to  the  defended  point.  Whenever  the 
turn  radius  is  less  than  the  distance  from  the  center  of 
curvature,  a  turn  is  possible  (Figure  3-3).  If  no  tum  is 

Kssible,  we  have  situation  three;  that  is,  no  bias  can 
used  (Figure  3-4).  Under  such  circumstances  a 
branch  to  the  data  storage  area  is  executed  and  stan¬ 
dard  values  for  prediction  are  used. 

If  a  turn  followed  by  a  linear  segment  is  to  be  used, 
the  process  continues.  Figure  3-5  presents  the  detailed 
geometry  of  the  tum.  The  new  heading  (0)  and  the 
time  to  turn  (T )  are  computed  using  the  following 
formulas: 


Rp  =  v/R02-  26  R0sin  H  (3.29) 


(R0  •  6  sin  H)  •  Rp  6  cos  H 

|  (3.30) 

Rp  (R0  -  6  sm  H)  ♦  6  *  cos  H 

dfl  *  A  ♦  H 

13.31) 

(3.32) 

Tt  •  0/6 

(3.33) 

6  m8g*  Ad 

(3.34) 

The  time  of  the  flight  to  the  end  point  of  the  turn  is 
compared  to  the  time  to  tum.  If  the  time  of  flight 
exceeds  the  time  to  (urn,  the  biased  prediction  is  used, 
otherwise  a  branch  to  the  storage  area  is  made  and  a 
standard  prediction  results. 

When  turn  bias  is  included,  values  for  the  velocity  in 
the  x  and  y  direction  are  computed.  The  prediction  in 
these  directions  is  a  linear  extrapolation  along  the  line 
from  the  end  of  the  turn  to  the  defended  point.  In  this 
case,  the  Z  coordinate  is  unbiased. 

For  all  cases,  the  last  action  is  a  branch  to  the  data 
storage  area  where  the  predicted  position  is  trans¬ 
formed  to  polar  coordinates  and  stored  along  with  the 
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Figure  3-2.  Straight-Line  Bias  Situation 


Figura  3-4.  No  Bias  Situation  (S  >  L) 


Figure  3-3.  Turn  to  Attack  Bias  Situation  (S  <  L) 


estimated  time  of  closest  approach  as  the  most  recent 
entry  in  the  PREDICTED  POSITION  File. 

Figure  3-6  summarizes  the  logical  structure  of  the 
Prediction  Event. 


3.4  GUN  MOTION  AND  FIRING  EVENT 


The  Function  of  this  event  is  to  compute  the  closest 
approach  of  a  shell  fired  at  current  clock  lime  given; 
the  history  of  predicted  positions,  the  ballistic  tables, 
and  the  actual  aircraft  flight  path.  Gun  servo  lag  and/ 
or  servo  regeneration  can  be  included.  When  included, 
lag  and/or  regeneration  characteristics  are  determined 
by  input  parameter  specification. 


In  the  following  discussion,  upper  case  letters  are 
used  to  designate  predicted  values,  subscripts  designate 
the  time  to  which  the  values  pertain,  primed  upper 
case  letters  are  used  to  indicate  actual  gun  pointing 
angles,  and  the  ietters  L  and  R  are  used  to  designate 
Lag  and  Regeneration  Standard  dot  notation  is  used  to 
indicate  rates  and  accelerations.  In  this  notation,  the 
azimuth  and  elevation  become: 


(3.35) 


Figure  3-5.  Geometry  of  a  Turn  Followed  by  a 
Linear  Segment  ■ 


The  lag  compulation  is  based  on  the  four  latest 
predicted  positions.  These  four  values  are  used  to 
estimate  current  rate  and  acceleration  values  for  the 
coordinate  in  question,  which  is  in  turn  used  in  the 
expression: 


L 


A“ 


Kv  Kg 


(3.37) 


d>' 


'*o  +  Lo  +  Ro 


(3.36) 


When  the  estimates  are  substituted  into  the  above 
expression  the  lag  equations  become: 


21 

1 

Ao  ♦ 

-13 

l»Mt 

2K.4,2 

2°*A 

-1 

A  x  . 

•17 

..  •>  1 

+  2K.V 

-At  + 

20KvAt  21CaAt2] 

|2QKyA|  +  2^^ 


A-2A, 

(338) 


Values  for  regeneration  are  computed  as  follows: 
The  current  predicted  position  is  transformed  to  carte¬ 
sian  coordinates.  These  are  combined  with  current 
values  for  smoothed  cartesian  rates  and  transformed 
back  to  polar  coordinates.  The  resulting  polar  rates  are 
used  to  estimated  polar  acceleration  using  the  follow¬ 
ing  equations: 


svV+yp2 

(3.39) 

S^XpXp  +  VpYpJ/S 

(3.40) 

b  -  -2.0  St) /D 

(3.41) 

•2.0  D  6$  -  S0Z 

0  r>2 

(3.42) 

direction.  The  magnitude  of  the  drop  vector  it  com¬ 
puted  as  -gt. 

Projectile  position  and  velocity  are  transformed  into 
cartesian  coordinates  for  the  computation  of  the  actual 
closest  approach  values. 

The  time  of  closest  approach  is  computed  by  correct¬ 
ing  the  estimated  time  of  closest  approach  which  is 
computed  by  the  predictor.  The  correction  is  made  by 
assuming:  that  the  aircraft  and  projectile  are  moving 
at  constant  rates  for  the  time  period  necessary  for  the 
projectile  to  get  from  either  the  predicted  position  to 
position  of  closest  approach,  or  in  the  case  when  the 
closest  approach  occurs  earlier  than  the  predicted  time, 
the  time  necessary  to  travel  from  the  closest  approach 
to  the  predicted  position.  Once  this  assumption  is 
made,  the  time  from  predicted  position  to  closest 
approach  can  be  computed  as: 

tc  =  |FxW-V  (3.45) 

where: 

T£  -  position  of  aircraft  relative  to  the  projectile  at 
predicted  impact. 

V  -  velocity  of  the  aircraft  relative  to  the  projec¬ 
tile  at  the  predicted  impact. 

The  time  of  closest  approach  is  then  given  by: 

t  =  tp  +  tc  (3.46) 

The  position  of  the  shell  at  closest  approach  is  given 
by: 

WVc  <3  47> 


These  values  are  used  in  the  standard  two-term  lag 
estimate  yielding: 


The  projectile  position,  at  current  time  plus  time  of 
flight  (estimated  time  of  closest  approach)  is  given  by: 
the  predicted  range  (D.).  the  adjusted  azimuth  (0',). 
and  the  adjusted  elevation  (^a).  The  projectile  velocity 
is  computed  as  the  vector  sum  of  the  remaining  veloc¬ 
ity  corr  -sp<  tiding  to  the  time  of  flight  and  a  drop 
compoi  ?nt  The  remaining  velocity  is  considered  to  be 
in  the  radial  direction  and  the  drop  in  the  vertical 


If  the  error  components  are  desired,  they  are  com¬ 
puted  as: 

E*K*Vtc  (3.48) 

With  R,  V,  and  t,  defined  as  above. 

The  time  of  closest  approach,  the  position  and  veloc¬ 
ity  of  the  shell  at  closest  approach,  the  position  and 
velocity  of  the  aircraft  relative  to  the  shell  at  closest 
approach,  and  the  radial  acceleration  of  the  aircraft  at 
closest  approach  arc  stored  in  the  SHELL  POSITION 
FILE. 

3.5  REPORT  GENERATOR  EVENT 

The  function  of  the  REPORT  GENERATOR 
EVENT  is  to  compute,  tabulate,  and  print  measures  of 
performance.  The  particular  performance  measures 
reported,  category  boundary  points  for  statistical  out¬ 
puts,  and  some  system  operating  parameters  are  input 
in  the  initialization  portion  of  this  event.  After  initial- 
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Figure  3-6  Logic  Flow  Diagram  of  the  Prediction 
Event 


ization,  the  event  occurs  once  for  each  replicate  of  an 
experiment.  On  this  entry,  requested  outputs  are  gener¬ 
ated,  printed  and  tabulated  for  an  across  replicate 
summary.  On  the  last  replicate  of  an  experiment,  the 
event  also  produces  a  report  of  across  replicate  statis¬ 
tics.  Figure  3-7  presents  the  logical  relationships  of  the 
above  functions. 

A  rather  unusual  feature  of  the  REPORT  GENER¬ 
ATOR  is  the  inclusion  of  the  affects  of  dispersion,  fire 
doctrine,  and  aircraft  vulnerable  volume  computations 
in  the  report  generating  loops.  This  arrangement  facili¬ 
tates  the  direct  comparison  (using  a  constant  aircraft 
path  and  stochastic  tracking  error  sequence)  of  differ¬ 
ent  dispersion  and  fire  doctrine  parameters.  The  pres¬ 
ence  of  a  ‘nest’  of  two  loops,  one  for  dispersion  and 
one  for  fire  doctrine,  allows  the  generation  of  reports 
for  up  to  IOC  combinations  of  dispersion  parameters 
with  fire  doctrine  parameters. 

In  the  following  sections,  the  content  and  method  of 
compulation  of  each  report  are  presented.  Reference  is 
made  to  the  following  utility  modules  which  perform 
major  support  functions: 

Aircraft  Position  Routine 

Fire  Doctrine  Routine  and 

Vulnerability  Routine 

3.5.1  Single  Replicate  Reports 

Measures  of  performance  generated  fall  into  four 
categories: 

a.  Miss  Distance  Measures. 

b.  Single  Shot  Probabilities. 

c.  Burst  Probabilities. 

d.  Course  Probabilities. 

Data  for  a  summary  report  is  always  generated  for 
each  category.  In  addition,  detailed  reports  may  be 
requested  in  any  combination  for  the  first  three 
categories. 

Figure  3-8  depicts  the  logical  structure  of  the  single 
replicate  portion  of  the  REPORT  GENERATOR.  Ta¬ 
ble  MI-2  presents  the  number  of  reports  which  can  be 
generated  for  one  replicate.  These  numbers  reflect  the 
loop  structure  of  the  REPORT  GENERATOR. 

3.5. 1.1  Miss  Distance  Reports 

Miss  distance  data  is  computed  using  values  stored 
in  the  SHELL  POSITION  FILE  and  accesses  to  the 
Aircraft  Position  Module.  Values  are  computed  with¬ 
out  considering  dispersion  or  fire  doctrine.  It  is  as¬ 
sumed  that  the  gun  fires  continuously  and  that  rounds 
fired  at  the  specified  point  in  time  have  zero  dispersion 
and  bias.  When  a  detailed  report  is  requested,  one  line 
of  output  for  each  sample  time  in  the  SHELL  POSI¬ 
TION  FILE  :s  generated.  Position  and  miss  distance 


3-8 


2037 1  186 

Figure  3-7.  Top  Level  Logic  of  the  Report  Generator 


values  are  for  the  point  of  closest  approach.  The 
following  data  is  presented: 

a.  Time  of  fire  (sec). 

b.  Time  of  flight  (sec). 

c.  Range  to  aircraft  (m). 

d.  Azimuth  of  aircraft  (degrees). 

e.  Elevation  of  aircraft  (degrees). 

f.  Radial  acceleration  of  aircraft  (g). 

g.  Range  error  (m). 

h.  Azimuth  error  (m). 

i.  Elevation  error  (m). 

j  Error  along  line  of  flight  (m).* 

k.  Error  perpendicular  to  line  of  flight  (m).* 

l.  RMS  error  (m).* 

NOTE:  *In  the  plane  normal  to  the  relative 

velocity  vector  at  closest  approach. 

A  summary  miss  distance  report  is  always  generated 
and  always  tabulated  for  reporting  in  the  across  repli¬ 
cate  report.  The  summary  report  consists  of  a  partition 
of  the  samples  computed  on  the  basis  of  RMS  error. 
Up  to  ten  partition  points  are  specified  and  the  number 


of  samples  which  do  not  exceed  the  partition  point  of 
each  category  is  tabulated  and  printed. 

Figures  3-9  and  3-10  present  samples  of  the  miss 
distance  report  formats. 

3.5. 1.2  8ingl«  Shot  Kill  Probability  Reports 

Single  shot  kill  probability  computations  are  based 
on  the  data  in  the  SHELL  POSITION  file,  the  disper¬ 
sion  parameter,  and  the  vulnerable  volume  of  the 
aircraft.  The  actual  computation  of  the  single  shot  kill 
probabilities  are  accomplished  in  the  vulnerability 
module.  The  module  is  invoked  once  for  each  entry  in 
the  SHELL  POSITION  FILE  (for  which  the  aircraft  is 
in  range;  and  returns  single  shot  hit  probabilities  for 
the:  fuselage  (P,).  the  wing  (P.).  and  the  fuselage-wing 
intersection  (P. ).  These  three  values  are  combined  using 
input  parameters  to  form  the  single  shot  kill  probabil¬ 
ity  for  the  aircraft  (P.).  The  most  usual  combination 
used  was: 

Pa  =  Pf+fVpi  •  <3-49> 

Values  of  P,  are  tabulated  for  summary  output  and  for 
the  across  replicate  report. 

When  a  detailed  report  is  requested,  one  line  of 
information  for  each  entry  in  the  SHELL  POSITION 
FILE  is  generated.  Each  line  contains  the  following 
information: 
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Figure  3-8.  Computations  For  a  Single  Replicate 


Tab)®  111*2.  Reports  Generated 


Report 

Number  of 

Reports 

Mils  Distance 

t 

Single  Shot  Kill  Probibility 

n 

Burst  Kill  Probability 

n  m 

Probability  Over  the  Course 

n  m 

n  *  number  of  dispersion  parameter  sets  s  10 
m  *  number  of  fire  doctrine  parameter  sets  s  10 

20871-534 


a.  Time  of  fire. 

b.  Probability  of  hitting  the  f uselage. 

c.  Probability  of  hitting  the  wings. 

d.  Probability  of  hitting  the  fuselage-wing  intersec¬ 
tion. 

e.  Single  shot  kill  probability. 

The  summary  report  consists  of  a  partition  of  the 
samples  on  the  basis  of  the  single  shot  kill  probability. 

Up  to  ten  partition  points  are  specified  in  the  input 
data.  The  number  of  samples  which  do  not  exceed  the 
partition  point  for  each  category  are  tabulated  and 
printed 

Figures  3-1 1  and  3-12  present  samples  of  the  single 
shot  output  reports. 

3.6. 1.3  Buret  Reports 

Burst  measures  are  computed  using  the  single  shot 
kill  probabilities  generated  in  the  single  shot  area  and 
the  fire  doctrine  parameters.  The  FIRE  DOCTRINE 
Module  is  accessed  to  determine  the  time  of  fire  of  the 
first  shot  of  a  burst,  the  duration  of  the  burst,  and  the 
number  of  rounds  (n)  in  the  burst.  This  data  is  used  to 
access  the  table  of  single  shot  kill  probabilities  and 
compute  an  average  value  over  the  burst  (p).  The 
probability  that  a  burst  kills  the  aircraft  (Pe)  is  com¬ 
puted  using  the  following  equation: 

PB  =  I  0  -  e  'nP  .  (3  50) 

When  requested,  a  detail  report  is  generated  The 
detail  report  contains  one  line  for  each  burst  which  can 

HISS  CIS 
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Figure  3-9.  Miss 


be  fired  according  to  the  fire  doctrine.  Each  line  con¬ 
tains  the  Following: 

a.  Time  of  fire. 

b.  Average  single  shot  kill  probability. 

c.  Probability  that  the  bunt  kills  the  aircraft. 

A  summary  burst  kill  report  is  always  generated  and 
always  tabulated  for  across  replicate  reporting.  The 
summary  consists  of  a  partition  of  the  sample  based  on 
the  computed  values  of  bunt  kill  probability.  The 
sample  is  partitioned  using  up  to  ten  input  category 
points. 

Figures  3-13  and  3-14  present  examples  of  the  burst 
kill  probability  reports. 

3. 6. 1.4  Over  the  Course  Meeeuree 

Two  measures  are  generated  over  the  entire  course. 
They  are  the  integral  of  the  single  shot  kill  probability 
over  time  which  is  independent  of  fire  doctrine  and  the 
kill  probability  over  the  course,  which  is  the  kill  proba¬ 
bility  resulting  from  combining  the  effects  of  all  bursts 
of  a  fire  doctrine. 

The  integral  of  the  single  shot  kill  probability  is 
computed  as  follows: 

j Psdt  *  X  At  (3  51) 


The  kill  probability  over  the  course  is  computed  as 
follows: 

Pc=l-e'nZPs  (3.52) 

where: 

n  «  number  of  rounds  per  burst 

P, "  average  single  shot  kill  probability  over  a 
burst,  and  the  sum  is  over  all  bursts  of  a  fire 
doctrine. 

3.5.2  Across  Replicate  Summary  Report 

The  across  replicate  summary  report  is  made  up  of 
the  summary  data  of  each  single  replicate  categories  of 
measures  and  the  kill  probability  over  the  course  aver¬ 
aged  over  the  set  of  replicates. 

Samples  of  these  reports  are  presented  in  Figures 
3-15  through  3-18. 
unce 
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Figure  3-10.  Miss  Distance  Detail  Report 
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Figure  3-11.  Single  Shot  Kill  Probebility  Detail  Report 
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Figure  3-12.  Single  Shot  Kill  Probability  Report 
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Figure  3-13.  Burst  Kill  Probability  Detail  Report 
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Figure  3-14.  Burst  Kill  Probability  Report 
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Figure  3-15.  Miss  Distance  Summary  Report 
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Figure  3*16.  Kill  Seconds  Over  the  Course  Report 
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Figure  3' 17.  Single  Shot  Kill  Probebility  Summery  Report 
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Figure  3-18.  Kill  Probebility  Over  the  Course  Report 
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SECTION  4 
UTILITY  PACKAGE 


4.1  AIRCRAFT  P08ITI0N  MODULE 

The  Aircraft  Position  Routine  provides  the  position 
and  velocity  of  the  aircraft  in  cartesian  coordinates 
centered  at  the  gun  for  the  simulation  time  specified  in 
the  invoking  command. 

Paths  are  represented  as  sequences  of  continuous 
path  segments  proceeding  from  an  initial  point  of 
known  position  and  velocity.  Each  path  segment  is 
either  a  segment  of  a  line,  a  segment  of  a  circle,  or  a 
segment  of  a  logarithmic  spiral.  All  path  segments  are 
in  the  horizontal  plane  coupled  with  motion  of  con* 
stant  acceleration  (possibly  zero)  in  the  vertical  direc¬ 
tion.  Smooth  positions  and  velocities  are  assured  by 
using  the  end  conditions  of  one  segment  is  the  initial 
conditions  of  its  successor.  The  end  point  of  a  segment 
is  determined  relative  to  its  initial  point  with  the 
‘duration’  of  the  segment  being  spcciAed  by:  length  for 
a  line  segment,  and  degrees  of  turn  for  both  circular 
and  spiral  segments.  The  only  position  and  velocity 
specification  made  is  for  the  initial  point.  Further  path 
description  is  in  the  form  of  the  following  statements: 

a.  Line  *X'  meters  in  length,  tangential  acceleration 
‘A,’  and  vertical  acceleration  *Af.’ 

b.  Circular  turn  ‘4’  degrees  right  or  left  radial 
acceleration  'A«'  and  vertical  acceleration  ‘A 

c.  Spiral  turn  ‘4’  degrees  right  or  left,  radial  accel¬ 
eration  'A*,'  tangential  acceleration  ‘At/  and 
vertical  acceleration  ‘A*/ 


For  linear  segments  position  and  velocity  are  deter¬ 
mined  according  to  the  equations: 


xo)-v(v0i  +  2  ATt2)  cos  e 

(4.1) 

Y(t)-Y0+(v0t  ♦  ^ATt2)sin0 

(4.2) 

)-Z0*V2ot  ♦  JaZ(2 

(4.3) 

V\(t)“  *(V0  4  Ajt)  coi  e 

(4.4) 

Vy(t)»(V0  4  Aft)  sin  9 

(4.5) 

Vz(t)-Vz  ♦  AZt 

(4.6) 

where:  6  ■  aircraft  heading 

For  circular  segments,  position  and  velocity  are 
determined  according  to  the  equations: 


R-V02/Ar 

(4.7) 

<I>* 

■ 

< 

50 

(4.8) 

«(t)-C0  +  $t 

(4.9) 

X(t)-Xo  ♦  R  [cos  90  -  cos  8(t)) 

(4.10) 

Y(t)  ■  Y0  ♦  R  (sin  60  4  ain  «(t)| 

(4.11) 

2(0 -Zo  *  VZ()t  4  *Azt2 

(4.12) 

V*(t)-  -V  sin  8(t) 

(4.13) 

Vy(t)  *  V  cos  9(t) 

(4.14) 

VZ(0  -  VZo  4  AZt 

(4.15) 

For  spiral  segments  position  and  velocity  are  deter* 
mined  according  to  the  equations: 


X*Aj/Ar  (4.16) 

V(t)-V0  ♦  ATt  (4.17) 

R-V02/AR  (4.18) 

*(t)-0o  -  Rloge  |v2(t)/2XV02j  (4.19) 


-  Icos  $0  -  cos  (0(t))l  J  (4.20) 


Y(t)-Y0 


V^t) 

Ar(  I  +4X^) 


cos  8(1)  -  coi  9 


o 


(4.21) 


4*1 


ZuWo  +VZot  +  jAztZ 

(4.22) 

Vjfltl  *  -V(i)  sin  I0(t)l 

(4.23) 

Vyjll*  V(t)cos  I0(t)) 

(4.24) 

VZ(t)-VZu  ♦  Azt 

(4.25) 

4.2  BALLISTICS  MODULE 

The  function  of  the  BALLISTICS  MODULE  it  to 
provide  shell  ballistics  information.  The  three  entry 
points  used  during  the  execution  phase  are  as  follows: 

a.  P RANGE.  Provides  range  as  a  function  of  time 
of  flight. 

b.  PVELOC.  Provides  remaining  velocity  as  a  func¬ 
tion  of  time  of  flight. 

c.  TRANGE.  Provides  time  of  flight  as  a  function 
of  range. 

During  the  initialization  phase,  a  table  of  remaining 
velocity,  range,  and  time  of  flight  is  input.  The  table 
may  be  ragged  and  need  not  start  at  the  desired  muzzle 
velocity.  Values  of  muzzle  velocity  and  a  ballistic 
coefficient  or  a  point  on  the  trajectory  are  also  input. 
The  input  tables  are  adjusted  to  fit  the  input  data  and 
interpolation  is  carried  out  to  produce  equally  spaced 
tables. 


tion  is  based  on  the  fact  that,  at  the  time  of  predicted 
impact,  the  range  to  the  aircraft  (R.)  and  the  range  to 
the  projectile  (R,)  must  be  equal.  In  cssense  then,  the 
Time  of  Flight  Module  numerically  solves  the 
equation: 

f(t)  ■  Ra(t)  -  Rp(t)  -  0  (4.26) 


for  the  smallest  value  of  t. 

The  values  of  aircraft  range  are  provided  by  the 
Prediction  Algorithm  Module  and  are  dependent  upon 
the  values  currently  stored  in  the  SENSED  DATA  rile 
as  well  as  the  particular  options  for  prediction  being 
simulated.  Values  for  R»(t)  are  obtained  by  invoking 
the  ‘range  as  a  function  of  time’  entry  of  the  ballistics 
package. 

A  start  value  for  the  solution  is  obtained  by  assum¬ 
ing  a  projectile  which  does  not  deaccelerate,  and  an 
aircraft  which  is  in  unaccelerated  motion: 


Tuan  u  **a(0) 


M) 


\/MO)*Vm2-V,2(0)j/(vm2  -  V,2(0)) 


where: 


(4.27) 


The  tables  produced  are  as  follows: 

a.  Range  as  a  function  of  Time  of  Flight. 

b.  Remaining  Velocity  as  a  function  of  Time  of 
Flight. 

c.  Time  of  Flight  as  a  function  of  Range. 

During  the  execution  phase,  any  values  required  are 
obtained  by  a  table  look-up  and  linear  interpolation  in 
the  appropriate  table.  Note  that  the  ballistic  informa¬ 
tion  is  isolated  in  this  module  and,  therefore,  any 
desired  ballistics  could  be  used  with  the  model.  Modi¬ 
fication  could  be  limited  to  performing  any  necessary 
initialization  in  the  initialization  entry,  thus  providing 
the  three  required  function:  by  using  the  appropriate 
names  and  calling  sequences. 

4  3  TIME  OF  FLIGHT  MODULI 

The  Time  of  Flight  Module  computes  the  lime  re¬ 
quired  for  a  projectile  to  travel  from  the  gun  to  a 
future  aircraft  position.  The  complexity  of  the  problem 
is  increased,  because  the  predicted  position  is  depen¬ 
dent  upon  the  lime  of  flight.  Closed-form  solutions  of 
the  problem  are  possible  only  under  some  rather  re¬ 
strictive  assumptions.  For  this  reason  numerical  meth¬ 
ods  have  been  employed.  The  lime  of  flight  computa¬ 


Vm  -  muzzle  velocity  of  gun. 

V,  •  current  aircraft  velocity. 

R.  -  current  range  rate  of  the  aircraft. 

An  initial  computation  of  the  function  f(t)  ■ 
R,(t)  •  R,(t)  is  made  using  the  start  value.  The  process 
is  continued  with  each  successive  computation  being 
made  at  a  new  time,  which  i*  the  current  value  incre¬ 
mented  by  a  function  of  the  start  value. 

Once  a  sign  change  has  been  obtained,  a  combina¬ 
tion  of  the  secant  method  and  interval  bisection  is 
employed  until  the  desired  accuracy  is  obtained.  In 
general,  no  more  than  five  computations  are  required 
to  derive  the  required  accuracy. 

4.4  PREDICTION  ALGORITHM  MODULE 

The  function  of  the  Prediction  Algorithm  Module  is 
to  extrapolate  current  available  aircraft  position  and 
velocity  data  to  a  specified  future  time  point,  (hereby, 
producing  an  estimate  of  aircraft  position  for  the 
specified  time.  The  module  has  the  capability  of  per¬ 
forming  the  extrapolation  by  a  number  of  methods, 
each  of  which  is  somewhat  flexible  in  that  its  specific 
behavior  is  determined  by  method  parameters.  The 
specific  method,  employed  during  a  simulation  run,  is 
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where: 


determined  by  the  designation  of  the  general  method 
;ype  and  its  required  parameters  during  the  initializa¬ 
tion  phase  of  the  run.  The  position  data  on  which  the 
module  operates  are  all  derived  from  the  output  of  the 
tensor  module.  These  data  arc  accessed  via  the  commu¬ 
nications  pool  (specifically  the  SENSED  DATA  File). 

Current  module  capabilities  include  the  following 
three  general  methods: 

a.  Quadratic  (includes  linear  by  parameter  specifi¬ 
cation). 

b.  Simple  polar. 

c.  Defense  of  a  known  point  (in  combination  with 
quadratic). 

4.4.1  Quadratlo  Prediction 

The  quadratic  prediction  it  made,  using  smoothed 
position  (X,  V.  Z).  velocity  (X,  Y,  Z).  and  acceleration 
(X.  Y.  Z)  values  computed  and  stored  in  the  SENSED 
POSITION  File  by  the  TRACK  Module.  The  predic¬ 
tion  it  computed  as  a  linear  estimate,  corrected  for 
accelerations  in  the  horizontal  plane,  and  changes  in 
total  aircraft  speed  induced  by  changes  in  aircraft 
altitude.  The  basic  estimate  for  a  time  increment  r  it 
given  by: 

X'-X+fcr  (4. 28) 

Y'mY  +  ir  (4.29) 

z  *  l  *  lT  (4.30) 

The  correction  for  acceleration  in  the  horizontal 
plane  it  made  whenever  the  absolute  value  of  the 
estimated  heading  rate  6  of  the  aircraft  exceeds  the 
value  of  an  input  threshold  parameter.  The  specifica¬ 
tion  of  an  excessively  large  threshold  value  leads  to 
suppression  of  the  acceleration  correction  and  linear 
prediction  results.  The  heading  rate  it  estimated  by  the 
equation: 

1  •  (S.Y  -  Xi)/(X2  *  X2)  <43,> 

When  the  correction  is  included,  the  resultant  pre¬ 
dictions  for  X  and  Y  are  at  follows: 


F  *  r2  var  ♦  b  (4.32) 


X'  •  X  +  Xr-  1/2  Ff  Y  (4,33) 


Y  ■  X*Yr-l/2FfX 


a  -  2.0  T„. 

b  -  2.0  (T„.  •  T,„  +  TP„J) 

Tp/,  -  Time  lag  of  position  due  to  constant 
velocity 

T,/,J  -  Time  lag  coefficient  of  position  due  to 
constant  acceleration 

T,„  -  Time  lag  of  velocity  due  to  constant 
acceleration 

The  correction  for  change  in  aircraft  speed,  due  to  a 
change  in  altitude,  is  made  whenever  the  absolute 
value  of  the  aircraft  velocity,  in  the  vertical  direction, 
exceeds  the  value  of  an  input  threshold  parameter,  as 
in  acceleration,  the  correction  can  be  suppressed  by 
specification  of  an  excessive  value  for  that  parameter. 
When  the  correction  is  included,  the  resultant  predic¬ 
tions  are  as  follows: 

X'  •  X  +  Xr-l/2  +  FldY  +  *(XZ/V2)l 

(4.35) 


V'  -  Y*Yr  +  l/2Fl£X-|(XZ/V2)] 

(4.36) 

V  -  Z  +  Zt  ♦  1/2  Fl«(Z 2/V2))  (4.37) 

where: 

F  is  defined  as  above 
g  is  acceleration  due  to  gravity 
V  is  the  aircraft  velocity. 

4.4.2  Simpi*  Polar  Prediction 

The  simple  polar  prediction  it  made  by  estimating 
the  range  rate  R,  the  azimuth  rate  A,  and  the  elevation 
rate  E.  Assuming  constant  rates,  and  using  linear 
extrapolation,  polar  coordinate  rates  are  estimated  by 
transforming  the  cartesian  coordinates  of  the  last  two 
entries  in  the  SENSED  POSITION  File  and  computing 
first-divided  differences.  Thus: 

i  •  (Bt  -  Bt-At)/A«  (4  38) 

A  *  (At  "  At-AtVAt  (4,39) 

i  •  (El  ~  Et-At m  (4.40) 


Estimated  range,  azimuth,  and  elevation  arc  given 
by  the  following  equations: 
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R  -  R,  +  Rr 

A  •  £,  +  A* 

F.  -  E,  ♦  Er 


(4.41) 

(4.42) 

(4.43) 


The  required  cartesian  estimates  are  obtained  by 
transforming  R',  A',  and  E'. 

4.4.3  D(fsns«  of  a  Known  Point 

The  majority  of  the  logic  and  computation  for  the 
defense  of  a  known  point  option  are  in  the  Prediction 
Event  Module.  However,  for  time  of  flight  computa¬ 
tions.  use  is  made  of  the  Prediction  Algorithm  Module. 
During  such  entries  to  the  module,  one  of  two  modes 
of  prediction  are  possible:  Quadratic  (discussed  previ¬ 
ously)  or  point-to-point  linear  interpolation.  An  indi¬ 
cator  switch,  set  in  the  Prediction  Event  Module,  con¬ 
trols  which  method  is  used  on  any  given  entry  to  the 
module.  The  following  pertains  to  point-to-point  linear 
interpolation,  used  for  a  portion  of  the  attack,  of  a 
known  point  option. 

Figure  4-1  presents  the  geometry  for  the  point-to- 
point  linear  interpolation.  The  point  Q  represents  the 
last  sensed  aircraft  position  and  P  the  point  being 
defended.  The  prediction  method  assumes  that  the 
aircraft  will  travel  along  the  line  connecting  the  points, 
at  rates  which  are  consistent  in  magnitude  to  the 
current  smoothed  aircraft  speed.  Distance  along  the 
line  is  computed  using  the  following  equations: 

AX  •  Xp  -  Xq  (4  44) 

AY  •  Yp  -  Y o  (4.45) 

AZ  ■  Zp  -  Zq  (4.46) 

D  •  /AX*"*  AY2  ♦  AZ*  (4.47) 


X  ^ 

lit)  ■  Vt  *  ^fl(t *  *  tT,  +  (4.48) 


Figure  4-1.  Geometric  Representation  of  the  Point- 
to-Point  Linear  Interpolation 

4.6  RANDOM  NUMBER  GENERATOR 

The  function  of  the  Random  Number  Generator  is 
to  produce  sample*  of  random  variables  from  either  a 
population  uniformly  distributed  on  the  interval  (0, 1  ], 
or  from  a  population  which  is  normally  distributed 
with  mean  zero  and  variance  one.  The  uniform  ran¬ 
dom  variables  are  generated  bv  using  an  additive 
congruous  process  with  31  seeds.  This  method  was 
based  on  work,  which  was  done  at  MIT’s  Lincoln 
Laboratory  in  1938,  by  B.F.  Oreen.  J.E.  Smith,  and  L. 
Klem.  The  sequence  of  samples  is  produced  by  using 
the  following  equations: 


Xj  •  (X).,  ♦  X,.n)Mod2'  J>n 
Rj-n*l  '  *j/2f 


(4.52) 


where: 


The  desired  estimate  is  then  obtained  by  interpola- 

n  - 

31 

lion: 

r  - 

29 

X  "  Xq  ♦  AX  Ut)/D 

(4.49) 

X,  - 

an  input  seed  variable 

Y  -  Yq  +  AY  U  D/D 

(4.30) 

X;  -  Xji  •> 

prestored  seed  values 

Z  •  Zq  ♦  AZ  Mt)/D 

(4,51) 

R,  - 

the  required  uniformly  distributed 
random  number*. 
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Normally-distributed  variables  are  computed  from 
the  uniformly-distributed  random  variables,  using  the 
well  known  Central  Limit  Theorem.  The  exact  equa¬ 
tion  employed  is  as  follows: 


where: 

N  -  the  required  sample 

R,  -  uniformly  distributed  random  variables 

This  module  has  the  ability  to  return  any  number  of 
normal  or  uniform  numbers  on  any  given  entry.  The 
sequence  of  samples  generated  is  a  function  of  the  first 
seed  value,  which  is  an  input  variable. 

4.6  COORDINATE  TRANSFORMATION 

There  are  two  coordinate  transformation  programs. 
One  converts  polar  coordinates  to  cartesian  coordinates 
and  the  other  converts  cartesian  coordinates  to  polar 
coordinates.  Both  programs  convert  the  three  position 
coordinates,  plus  their  respective  rates.  The  transfor¬ 
mations  are  not  standard,  due  to  the  fact  that  the 
azimuth  angle  is  measured  counter-clockwise  from  the 
positive  y-axis.  This  system  was  used  in  an  attempt  to 
remain  as  compatible  as  possible  with  systems  em¬ 
ployed  in  the  University  of  Michigan  Gun  Model. 

The  cartesian  to  polar  transformation  uses  Equations 
4.5b  through  4.61.  The  equations  are  computed  in  the 
order  shown  and  earlier  results  are  used  in  later  com¬ 
putation  in  order  to  minimize  computer  time. 


R  •  v/x2  t  Y2  +  Z2 

(4.55) 

S  ■  \/x2  +  Y2 

(4,561 

6  -  Tan*'  (-X/Y) 

(4.57) 

0  ■  Tan'*  (Z/S) 

(4.58) 

ft  •  (XX  ♦  YY  +  ZZ)/R 

(4.59) 

$  -  (XY  -  XY)/S2 

(4.60) 

0  •  (Z  -  Zft/R)/S 

(4.61) 

The  polar  to  cartesian  conversion  uses  the  following 
equations: 


S  •  R  cos  0 

(4.62) 

Z  ■  R  sin  0 

(4.63) 

X  ■  -S  sin  9 

(4.64) 

Y»Scos  9 

(4.65) 

T  -  ft/R  •  Z0/S 

(4.66) 

*»XTYrf 

(4.67) 

Y  •  X6  +  YT 

(468) 

Zft/R  +  $0 

(469) 

4.7  FIRE  DOCTRINE  MODULE 

The  function  of  the  Fire  Doctrine  Module  is  to 
provide  a  FIRE  -  NO  FIRE  decision  for  a  given  point 
in  time.  When  a  FIRE  decision  is  made,  Fire  Doctrine 
Module  provides  values  for.  burst  duration,  rate  of  fire, 
and  number  of  rounds  in  the  burst. 

The  current  configuration  works  in  conjunction  with 
the  report  generator,  which  it  responsible  for  main¬ 
taining  the  count  of  rounds  available  for  firing  and  not 
initiating  the  Fire  Doctrine  Module,  unless  there  are 
rounds  on  the  mount.  The  following  parameters  are 
input  during  initialization: 

a.  Open  fire  lime.  First  time  at  which  a  burst  may  be 
fired. 

b.  Round  per  burst.  Number  of  rounds  in  one  burst. 

c.  Burst  duration.  Period  of  time  over  which  burst 
will  be  fired. 

d.  Dead  time.  Minimum  time  between  firing  bursts. 

e.  Load  time.  Time  required  to  reload  the  mount. 

f.  Total  rounds.  Number  of  rounds  on  the  mount. 

The  Fire  Doctrine  Module  logic  and  REPORT 
GENERATOR  EVENT  logic  work  'ogether  to  assure 
that  bursts  are  fired  in  a  manner  consistent  with  the 
input  parameters.  If  the  number  of  rounds  in  the 
mount  is  greater  than  the  rounds  in  the  burst,  the  logic 
spaces  the  bursts  by  the  dead  time.  When  ammunition 
is  exhausted,  the  logic  prevents  the  gun  from  firing, 
until  sufficient  time  has  elapsed,  to  accomplish  reload¬ 
ing  the  mount. 

4.8  VULNERABILITY  MODULE 

The  function  of  the  Vulnerability  Module  is  to  oper¬ 
ate  on  the  data  in  the  SHELL  POSITION  FILE  and 
produce  values  for  the  single  shot  hit  probability  for 
the  fuselage,  wing,  and  fuselage-wing  intersection. 

The  equation  used  to  compute  the  single  shot  hit 
probability  is  as  follows: 
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(4.70) 


where: 


|  jjl/2  - 1/2  Mj  W+S)  ’lM, 
Pss  |j  +  S|,/2' 


where: 

(JJ  is  ihe  inverse  of  the  matrix  JV], 

|V]  describing  the  diffuse  ellipsoid  represent¬ 
ing  the  aircraft  segment. 

[S]  is  the  convariance  matrix  of  the  round  to 
round  dispersion  of  the  shots. 

M,  is  the  vector  of  displacement  of  the  aim 
point  (point  of  closest  approach)  from  the 
center  of  the  target. 

The  following  paragraphs  expand  the  above  defini¬ 
tion  in  terms  of  aircraft  and  projectile  parameters. 

4.8.1  The  J  Matrix 

The  J  matrix  is  the  inverse  of  the  V  matrix,  which 
represents  the  projection  of  the  ellipsoid,  representing 
the  fuselage,  the  wing,  or  the  fuselage-wing  intersec¬ 
tion  onto  a  plane  normal  to  the  relative  velocity  vector 
of  the  shell  with  respect  to  the  aircraft.  The  V  matrix  is 
obtained  by:  1)  rotating  the  aircraft  (actually  the 
appropriate  matrix  representing  the  ellipsoid  of  vul¬ 
nerable  volume  corresponding  to  the  aircraft  compo¬ 
nent  in  question)  to  an  orientation  in  space  corre¬ 
sponding  to  the  true  position,  heading,  and  bank  angle 
of  (he  aircraft  at  the  time  in  question:  2)  project  this 
rotated  ellipsoid  onto  the  plan  normal  to  the  relative 
velocity  vector  of  the  shell  with  respect  to  the  aircraft 
at  the  point  where  the  ellipse  is  of  maximum  extent  in 
the  U  direction  to  obtain  the  required  V  matrix.  The 
above  process  in  equation  form  is  as  follows: 


Define: 


Ml) 

M21 

m3) 

M  * 

m21 

M22 

M32 

M3I 

M32 

m33 

(4.71) 


H 


-2  0 


0  p  0 

0  0-2 
c 


is  the  aircraft  ellipsoid 


(4.72) 


R  »  BT^TAR  |eR  +  tt/2J  (4.73) 


B 

e 

<t> 

A* 

e„ 

Then: 


V 


Rotation  Matrix  for  Bank  angle. 
Rotation  Matrix  for  Climb  angle. 
Rotation  Matrix  for  Heading  angle. 
Rotation  Matrix  for  Relative  azimuth. 
Rotation  Matrix  for  Relative  elevation. 


M 


M31M31 


11  '  ~M 


33 


M21 


m32m3, 


M 


33 


m21 

m22 


M31M32 


M32M32 

M33 


J  =  V1 

4.8.2  The  Covariance  Matrix  8 


(4.74) 


The  covariance  matrix  [S)  is  computed  as  the  sum  of 
the  matrix  | <r,2)  of  angular  dispersion  and  the  matrix 
[<rj  J  or  muzzle  velocity  dispersion. 

The  angular  dispersion  matrix  (<r,2)  is  obtained  by 
rotation  of  the  angular  dispersion  matrix,  in  gun  coor¬ 
dinates  by  an  e'evation  angle,  e,  and  an  azimuth  angle, 
A.  3nd  then  back  through  a  relative  azimuth  angle,  A,, 
and  an  augmented  elevation  angle  (e,  +  n/2).  Only 
the  upper  right  hand  four  elements  arc  retained  in 
equation  form: 
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[o.2]  •  c[»,2]  cT  <4”> 

CT  »  CTAAf[er  +  w/2)  (4.76) 


The  muz2le  velocity  dispersion  matrix  Jo-*,*]  is  com¬ 
puted  by  treating  the  error  as  an  error  along  the  line  of 
flight,  due  to  an  early  or  late  arrival  of  the  shell,  and 
applying  a  rotation  to  obtain  the  appropriate  values  in 
the  plane  normal  to  the  shell  trajectory  near  the  air¬ 
craft.  The  error  is  estimated  using  a  3/2  law  for 
ballistics  and  the  rotation,  R,  is  the  same  as  R  used  in 
the  definition  of  the  J  Matrix. 

In  equation  form: 


where: 

o\o  m  muzzle  velocity  dispersion 
V0  -  muzzle  velocity 

V  «•  aircraft  velocity 

V„  -  projectile  velocity 

De  -  range  of  projectile 


0 


1° 


0 


0 


(4.78) 

4.8.3  Tha  M,  V actor 

The  M,  vector  is  the  displacement  from  the  point  of 
closest  approach  to  the  target  center.  The  vector  is 
composed  of  the  error  due  to  prediction  and  the  error 


due  to  any  systematic  difference  between  the  muzzle 
velocity  used  for  prediction  and  the  actual  muzzle 
velocity- 

M,  »  Mp  +  Mb  (4.79) 

where: 

Mp  ■  displacement  due  to  prediction  error 
Mb  -  displacement  due  to  muzzle  velocity  bias. 

Mp  is  obtained  by  rotation  of  the  prediction  errors 
from  the  gun  based  system  to  the  projectile  based 
system  and  retaining  only  the  U  and  V  values. 

Mp  =  [XYZ]  [Ar]  [Ef  ♦  7T/2)  (4.80) 

Mb  is  produced  by  computing  the  error  in  gun 
coordinates  and  rotating  it  to  projectile  coordinates. 
The  error  computation  is  an  error  due  to  early  or  late 
arrival  assuming  a  3/2  ballistics  law  while  the  projec¬ 
tile  is  in  the  neighborhood  of  the  target.  In  equation 
form: 


Mg  =  [0  Bv  0)  R  (4.82) 


where: 

V,  ■  aircraft  velocity 
Vp  -  projectile  velocity 
AV,  -  muzzle  velocity  bias 
Vp  -  muzzle  velocity 
Dp  ■  range  of  the  projectile 

R  »  the  rotation  matrix  defined  in  Section 

3.8.1. 
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